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Given  a  physics-based  problem,  we  often  have  choices  in  deciding 

a)  which  Reduced-Order  Model  to  use 

b)  which  mathematical  method  to  use  to  analyze  the  ROM 


Examples  come  from  Ocean  Modeling,  Electromagnetism,  and 
Structural  Mechanics 


Typical  Problem  Genesis: 


•  Adhesively  Bonded  Joints 


Adhesive 


Adherends 


Aerospace  Sandwich  Structures 


Boundary  Value  Problem: 

Mixed  form  in  terms  of  Stress,  Infinitesimal 

Strain,  and  Displacement 


Boundary  Value  Problem: 
Expressed  in  terms  of  Displacement 


Boundary  Value  Problem: 
Expressed  in  terms  of  Displacement 

Potential 


Boundary  Value  Problem: 
Expressed  in  terms  of  Airy's  Stress  Function 

and  Solution  Ansatz 


d>(jc,z)  =  A0{x)  +  zAl(x)  +  ^Fn(*)sin| 

n=l 


Example  Loading  Case:  Displacements, 
Geometry,  and  Properties 


Spectral-Collocation  Analysis  Results: 
Mid-plane  Shear  Stress  oxz(x,n/2)  and 
Interfacial  Shear  Stresses  oxz(x,0) 


COMSOL  Structural  Mechanics  Mesh: 
Shown  in  vicinity  of  stress  free  surface 


COMSOL  Structural  Mechanics 
2D  Plane-Stress  Analysis  Results: 
Mid-plane  Shear  Stress  axz(x,r|/2) 


sxy  shear  stress  global  sys.  [MPa] 


x 


COMSOL  Structural  Mechanics 
2D  Plane-Stress  Analysis  Results: 
Interfacial  Shear  Stress  cjxz(x,0) 


sxy  shear  stress  global  sys.  [MPa] 


x 


Back  to  Displacement  BVP 


VT  =  F 


Use  COMSOL’s  General  PDE 
Solver 


Gamma  is  a  2  x  2  tensor 


COMSOL  PDE  (General  Form)  Solver  Mesh: 
Shown  in  vicinity  of  stress  free  surface 


sky  shear  stress  [MPa) 


COMSOL  PDE  Solver  Analysis  Results: 
Mid-plane  Shear  Stress  axz(x,r|/2) 


sky  shear  stress  [MPa) 


0  - 1 - 1 - 1 - 1 - 1 - 1 - 1 - L 

012345678 


*[m]  xicr3 


sky  shear  stress  [MPa) 


COMSOL  PDE  Solver  Analysis  Results 
Interfacial  Shear  Stress  cjxz(x,0) 


sky  shear  stress  [MIFa| 

0.2  - 1 - 1 - 1 - 1 - r 


0.L8 


0.16 


0„L4 
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012345678 


Some  Conclusions/Observations 


The  Spectral  Collocation  method  describes  a  shear  stress  that  does  not  have 
a  singularity  at  the  corners;  this  seems  to  be  the  expected  result  from  a 
mechanical  point  of  view. 


The  Structural  Mechanics  result  seems  to  suffer  from  artificially  large  singularities 
at  the  corners;  did  we  implement  it  poorly? 


The  General  PDE  solver  seems  to  be  a  natural  way  to  pose  this  pose  this  problem 
In  COMSOL.  It  gives  good  result,  although  it  seems  to  still  suffer  from  some  level  of 
Numerical  singularity  at  the  corners. 


